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Abstract. We consider the zeros on the boundary dfl of a Neumann eigenfunction (p\ of 
a real analytic plane domain f2. We prove that the number of its boundary zeros is 0(A) 
where —Aipx = A^(/?a- We also prove that the number of boundary critical points of either a 
Neumann or Dirichlet eigenfunction is 0(A). It follows that the number of nodal lines of (p\ 
(components of the nodal set) which touch the boundary is of order A. This upper bound is 
of the same order of magnitude as the length of the total nodal line, but is the square root 
of the Courant bound on the number of nodal components in the interior. More generally, 
the results are proved for piecewise analytic domains. 



This article is concerned with the high energy asymptotics of nodal lines of Neumann 
(resp. Dirichlet) eigenfunctions ipx on piecewise real analytic plane domains Q C M^: 



Here, is the interior unit normal. We denote by {fx^} an orthonormal basis of eigen- 
functions of the boundary value problem corresponding to eigenvalues Aq < Ai < A2 ■ • • 
enumerated according to multiplicity. The nodal set 



is a curve (possibly with self-intersections at the singular points) which intersects the bound- 
ary in the set N"^^ fl dQ of boundary nodal points. The motivating problem of this article is 

the following: how many nodal lines (i.e. components of the nodal set) touch the boundary? 
Since the boundary lies in the nodal set for Dirichlet boundary conditions, we remove it from 
the nodal set before counting components. Henceforth, the number of components of the 
nodal set in the Dirichlet case means the number of components of N'^^\dQ. 

Theorem 1. Let fl be a piecewise analytic domain and let nQQ{Xj) be the number of com- 
ponents of the nodal set of the jth Neumann or Dirichlet eigenfunction which intersect dQ. 
Then ngn{Xj) = 0{Xj). 

For generic piecewise analytic plane domains, zero is a regular value of all the eigenfunc- 
tions (fx, i.e. Vipx 7^ on J\f^^ [U]; we then call the nodal set regular. Each regular nodal set 
decomposes into a disjoint union of connected components which are homeomorphic either 
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1. Introduction 




(1.1) 



{x eQ: (px{x) = 0} 
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to circles contained in the interior ^1° of or to intervals intersecting the boundary in two 
points. We term the former 'closed nodal loops' and the latter 'open nodal lines'. Thus, we 
are counting open nodal lines. Such open nodal lines might 'percolate' in the sense of [Ze] 
(i.e. become infinitely long in the scaling limit Q XQ), or they might form -'small' 
half-loops at the boundary. Our methods may be useful in counting each type of component 
and it seems an interesting direction for future work. 

For the Neumann problem, the boundary nodal points are the same as the zeros of the 
boundary values ipx\dn of the eigenfunctions. The number of open nodal lines is thus twice 
the number of boundary nodal points. Hence we can count open nodal lines by counting 
boundary nodal points. In the Neumann case, our result follows from: 

Theorem 2. Suppose that Q C M."^ is a piecewise real analytic plane domain. Then the 
number n(Xj) = i^J^ip^. H dQ of zeros of the boundary values ip\-\dn of the jth Neumann 
eigenfunction satisfies nlXj) < CXj, where C is a constant depending only on Q. 

This is a more precise version of Theorem [1] in cases such as integrable billiard domains 
(rectangles, discs, ellipses) where the entire nodal set is connected due to the large grid of 
self- intersection points of the nodal set. The analogous result in the Dirichlet case is stated 
in Corollary |H Counting boundary nodal points of eigenfunctions has obvious similarities 
to measuring the length of the interior nodal line, and our results show that the order of 
magnitude is the same. We recall that S. T. Yau conjectured that in all dimensions, the 
hypersurface volume should satisfy cXj < 7Y"~^(A/'^J < CXj for some positive constants 
c, C depending only on (M, g) \Y1\ IY2j . The lower bound was proved in dimension two 
for smooth domains by Briining-Gomes [BG] and both the upper and lower bounds were 
proved in all dimensions for analytic {M,g) by Donnelly- Fefferman |DFl IDF2] (see also [L]). 
For general C°° Riemannian manifolds (M, g) in dimensions > 3 there are at present only 
exponential bounds [HHLl IHS] . Our methods involve analytic continuation to the complex 
as in [DFt IDF2t [L] , and it is not clear how to extend them to C°° domains. 

In comparison to the order 0{Xj) of the number of boundary nodal points, the total 
number of connected components of A/"^^ has the upper bound 0(Aj) by the Courant nodal 
domain theorem. Only in very rare cases is it known whether this upper bound is achieved 
(in terms of order of magnitude). When the upper bound is achieved, the number of open 
nodal lines in dimension 2 is of one lower order in Xj than the number of closed nodal 
loops. This effect is known from numerical experiments of eigenfunctions and random waves 
|BGS[ IFGSj . The only rigorous result we know is the recent proof in |NSj that the average 
number of nodal components of a random spherical harmonic is of order of magnitude A^. In 
special cases, the number of connected components can be much smaller than the Courant 
bound, e.g. two or three for arbitrarily high eigenvalues [Lew]. 

Our methods also yield estimates on the number of critical points of (px- which occur on 
the boundary. We denote the boundary critical set by 



rivative automatically equals zero on the boundary, while in the Dirichlet case q E C^ 
dyipx.{q) = since the boundary is a level set. 



{qedn:{dvx,){q) = 0}. 



In the case of Neumann eigenfunctions, q E 
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A direct parallel to Yau's conjecture for interior critical points of generic analytic metrics 
would be the Bezout bound ij^C^^^ < CA" in dimension n, and A"~^ for boundary critical 
points. However, this bound is unstable since the critical point sets do not even have to 
be discrete when the eigenfunctions have degenerate critical points, and the true count in 
the discrete case might reflect the size of the determinant of the Hessian at the critical 
point. Note that there is no non-trivial lower bound on the number of interior critical points 
[JNj . Related complications occur for boundary critical points. For instance, the radial 
eigenfunctions on the disc are constant on the boundary; thus, boundary critical point sets 
need not be isolated. We therefore need to add a non-degeneracy condition on the derivative 
dt{^Px\dn) to ensure that its zeros are isolated and can be counted by our methods. 

We phrase the condition in terms of the Pompeiu problem and Schiffer conjecture, which 
asserts the disc is the only smooth plane domain possessing a Neumann eigenfunction which 
is constant on the boundary. In pBerj) it is proved that the disc is the only bounded simply 
connected plane domain possessing an infinite sequence of such Neumann eigenfunctions. We 
say that the Neumann problem for a bounded domain has the asymptotic Schiffer property 
if there exists C > such that, for all Neumann eigenfunctions (p\, 

WdtfxWman) ^ ^_cA 2) 
||</'A||L2(an) ~ 

Here, the norms refer to the restrictions of the eigenfunction to dQ. It seems plausible 
that Berenstein's result might extend to this asymptotic Schiffer property, i.e. that the disc 
is the only example where fll.2p fails for an infinite sequence. 

Theorem 3. Let Q cM."^ be piecewise real analytic. Suppose that v^a Ian satisfies the asymp- 
totic Schiffer condition U.S\) in the Neumann case. Then the number of n ^^.j^^i^Xj) = i^C^p^^ 
of critical points of a Neumann or Dirichlet eigenfunction ipx. which lie on dfl satisfies 
^crit^^j) — /^'^ some C depending only on Q. 

This is apparently the first general result on the asymptotic number of critical points 
of eigenfunctions as Xj —>■ oo. Our results on boundary critical points give some positive 
evidence for the Bezout upper bound in the analytic case, and simple examples such as 
the disc show that there does not exist a non-trivial lower bound (see §2]). In general, 
counting critical points is subtler than counting zeros or singular points (i.e. points where 
^x,{x) = d^x,{x) = 0; see [HSIEHL]). 

In the case of Dirichlet eigenfunctions, endpoints of open nodal lines are always boundary 
critical points, since they must be singular points of v^a^ . Hence, an upper bound for ^crit('^i) 
also gives an upper bound for the number of open nodal lines. 

Corollary 4. Suppose that C zs a piecewise real analytic plane domain. Let nQ^{Xj) 
be the number of open nodal lines of the jth Dirichlet eigenfunction, i.e. connected compo- 
nents of {(fx^ = 0} C Q° whose closure intersects dfl. Then nan(Aj) = 0{Xj). 

The question may arise why we are so concerned with piecewise analytic domains C M^. 
By this, we mean a compact domain with piecewise analytic boundary, i.e. dQ is a union of 
a finite number of piecewise analytic curves which intersect only at their common endpoints 
(cf. |HZ] ). Our interest in such domains is due to the fact that many important types of 
domains in classical and quantum billiards, such as the Bunimovich stadium or Sinai billiard. 
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are only piecewise analytic. Their nodal sets have been the subject of a number of numerical 
studies (e.g. |BGSt IFGS] ). In |TZ] we consider piecewise analytic Euclidean plane domains 
with ergodic billiards, which can never be fully analytic. 

The results stated above are corollaries of one basic result concerning the complex zeros and 
critical points of analytic continuations of Cauchy data of eigenfunctions. When dQ G C^, 
the eigenfunctions can be holomorphically continued to an open tube domain in projecting 
over an open neighborhood W in of Q which is independent of the eigenvalue. We 
denote by flc C the points ( = x + iC, G with x G Then ipx^ (x) extends to a 
holomorphic function </3^^(C) where x E W and where |^| < eg for some eo > 0. We mainly 
use the complexifications to obtain upper bounds on real zeros, so are not concerned with 
the maximal eo, i.e. the 'radius of the Grauert tube' around dQ, and do not include the 
radius in our notation for the complexification. 

Assuming dQ real analytic, we define the (interior) complex nodal set by 

and the (interior) complex critical point set by 

C^^^={CeQc:d^l{O = 0}. 

We are mainly interested in the restriction of (p'^. to the complexification {dQ)c of the 
boundary, i.e. the open complex curve in obtained by analytically continuing a real 
analytic parameterization Q : ^ dfl. The map Q admits a holomorphic extension to an 
annulus A{e) (see (13.11) ) around the parameterizing circle and its image Qc{A{e)) C 
is an annulus in the complexification of the boundary; it is analogous to a Grauert tube 
around the real analytic boundary in the sense of [GSll ILSl] . We then define the boundary 
complex nodal set by 

and the (boundary) complex critical point set by 

Cj;;^- = {CG9fic:rf<.(C) = 0}. 

More generally, we may assume dfl is piecewise real analytic and holomorphically extend 
eigenfunctions to the analytic continuations of the real analytic boundary arcs. The radii of 
these analytic continuations of course shrink to zero at the corners. 

Theorem 5. Suppose that Q C M."^ is a piecewise real analytic plane domain, and denote by 
{dQ)c the union of the complexifications of its real analytic boundary components. 

(1) Let n{Xj,dQc) = ij^M^^^. Then, n{Xj,dQc) = 0{Xj). The 0-symbol is uniform in 
the radius of {dQ)c- 

(2) Suppose that the Neumann eigenfunctions satisfy and let n (.^.j^f^Xj , dQc) = 
i^C^^c, Then, n,^,i(A„9fic) = 0(A,). 

The theorems on real nodal lines and critical points follow from the fact that real zeros 
and critical points are also complex zeros and critical points, hence 



n{Xj) < n{Xj, dQc); ^crit(^i) ^ ^crit(^i' ^^c)- 



(1.3) 
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All of the results are sharp, and are already obtained for certain sequences of eigenf unctions 
on a disc (see §2]). If the condition (11. 2p is not satisfied, the boundary value of Lp\ is o(e~'"'^) 
close to a constant for all A. It is very likely that this forces the boundary values to be 
constant, but it would take us too far afield in this article to prove it. 

Although our main interest is in counting open nodal lines, the method of proof of Theorem 
|5] generalizes from dVt to a rather large class of real analytic curves C C f2, even when dVt is 
not real analytic. Let us call a real analytic curve C a good curve if there exists a constant 
a > so that 

\\V\\\L-i(C) 

Here, the norms refer to the restrictions of the eigenfunction to C and to dVt. The 
following result deals with the case where C C dVL is an interior real-analytic curve. The 
real curve C may then be holomorphically continued to a complex curve C obtained 
by analytically continuing a real analytic parametrization of C. 

Theorem 6. Suppose that dM? is a C°° plane domain, and let C d Vt he a good interior 
real analytic curve in the sense of \l-4\) - Let n{Xj,C) = ^A/"^^ (1 C be the number of 

intersection points of the nodal set of the j-th Neumann (or Dirichlet) eigenfunction with C . 
Thenn{Xj,C) =0{Xj). 

Although the upper bounds are sharp for some domains, we do not present necessary or 
sufficient conditions on a domain that the bounds on zeros or critical points are achieved on 
that domain for some sequence of eigenfunctions. We do not know any domain for which 
they are not achieved, but there are few domains where the bounds can be explicitly tested. 
The boundary (or rather its unit ball bundle) is naturally viewed as a kind of quantum 'cross 
section' of the wave group [HZj . The growth rate of the modulus and zeros of Cauchy data 
of complexified eigenfunctions depend on what kind of 'cross section' the boundary provides. 
In work in progress |TZj . we show that at least for some piecewise analytic domains with 
ergodic billiards, the the number of complex zeros of '^'x^ldUc is ~ CXj. It seems that this 
asymptotic refiects the fact that the boundary is a representative cross section in this case. 

We note that some of the methods and results of this paper are restricted to dimension 
two. In higher dimensions, zeros of the Cauchy data are not isolated and we would have to 
count numbers of components of the boundary nodal set. This seems inaccessible at present. 

The organization of this article is as follow: In ^ we use the layer potential representations 
of Cauchy data of eigenfunctions, or equivalently the representation in terms of the Calderon 
projector, to analytically continue eigenfunctions. The analytic continuation of the layer 
potential representation has previously been studied by Vekua [Vj, Garabedian [G], and in 
the form we need by Millar \M1\ \M2\ IM2j . The analytic continuation is somewhat subtle 
due to the presence of logarithms in the layer potentials, and does not appear to be well- 
known; so we present complete details (which are sometimes sketchy in the original articles). 
In ^ we relate growth of complex zeros to growth of the log modulus of the complexified 
eigenfunctions. In ^ - ^ we prove the main results. The complexified layer potential 
representation is used to obtain an upper bound on the growth rate of the complexified 
eigenfunctions y?^ in a fixed complex tube around the boundary as A — > cxd. The estimate 
is simpler for interior curves (^ since on the boundary the analytic continuation involves a 
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Volterra operator that must be inverted. Almost the same method gives analogous results 
on critical points; for the sake of brevity, the argument is only sketched in ^ 

We would like to thank L. Ehrenpreis, C. Epstein, P. Koosis and M. Zworski for informative 
discussions of analytic continuations of eigenfunctions and logarithmic integrals. 

2. Examples 

We begin with examples illustrating some of the issues we face. Eigenfunctions are only 
computable in (quantum) completely integrable cases, and at present the only known ex- 
amples are the unit disc, ellipses and rectangles. It is a classical conjecture of Birhoff that 
ellipses are the only smooth Euclidean plane domains with integrable billiards, so one does 
not expect further explicitly computable examples. In addition, one can construct approxi- 
mate eigenfunctions, or quasi-modes for many further domains |BB] ; it is plausible, although 
it is not proved here, that our results extend to real analytic quasi-modes. 

2.1. The unit disc D. The standard orthonormal basis of real valued Neumann eigefunc- 
tions is given in polar coordinates by y9m „(r, = Cm,,n sin m6'Jm(j^ ^r), (resp. C m,n cos m9Jm{j'„ 
where j'^^ is the nth critical point of the Bessel function and where Cm,n is the nor- 
malizing constant. The eigenvalue is A^^ = {j'mnf'- The parameter m is referred to as the 
angular momentum. Dirichlet eigenfunctions have a similar form with j'^ ^ replaced by the 
nth zero im,n of Jm- Nodal loops correspond to zeros of the radial factor while open nodal 
lines correspond to zeros of the angular factor. 

If we fix m and let \m,n oo we obtain a sequence of eigenfunctions of bounded angular 
momentum but high energy. In the sin case (e.g.), the open nodal lines consist of the union 
of rays Cm = {6 = k = 0, . . . , -^^^j through the mth roots of unity. Hence, for each 
m there exist sequences of eigenfunctions with A ^ oo but with m open nodal lines; hence, 
there exists no lower bound on the number of nodal lines touching the boundary in terms of 
the energy. This example also shows that there cannot exist a general unconditional result 
counting intersections of nodal lines with interior curves, since (pm,n\cm = hence the 
'number' of nodal points on the interior curve Cm is infinite. In particular. Cm is not 'good' 
in the sense of (11.41) . 

At the opposite extreme are the whispering gallery modes which concentrate along the 
boundary. These are eigenfunctions of maximal angular momentum (with given energy), 
and Am ~ fn. As discussed in |BB] . they are asymptotically given by the real and imaginary 
parts of e"^^"^^ Aip{p~^^^ Xm^ y) . Here, Aip{y) := Ai{—tp + y) where Ai is the Airy function 
and {—tp} are its negative zeros. Also, s is arc-length along dD, p is a normalizing constant 
and y = 1 — r. Whispering gallery modes saturate the upper bound on the number of open 
nodal lines. 

2.2. An ellipse. We express an ellipse in the form -|- = 1, < a < 1, with foci at 
{x, y) = (±a, 0). We define elliptical coordinates (</?, p) by (x, y) = (a cos (p cosh p, a sin (p sinh p). 
Here, < p < pmax = cosh~^a~^, < cp < 2tt. The lines p = const are confocal ellipses 
and the lines = const are confocal hyperbolae. The foci occur at ^9 = 0, vr while the origin 
occurs at p = 0, = f . 
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The eigenvalue problem separates into a pair of Mathieu equations, 

r dlG - cos'^ = -\^G 

{ (2.1) 
[ d^F - c2 cosh^ pF = X^F 

where c is a certain parameter. The eigenfunctions have the form \E'm,n(v', p) = Cm,nFm,n{p) ■ 
Gm,n{'^) where, Fm,n{p) = Gcmip, ^) and Gm,n{^) = ccmi^p, ^) (and their sin analogues). 
Here, cCm, Gcm are special Mathieu functions (cf. [C] (3.10)-(3.2)). The Neumann or Dirich- 
let boundary conditions determine the eigenvalue parameteres fc„c. The nodal lines are of 
course given by {G = 0} U {F = 0}. For more details and computer graphics of elliptic 
bouncing ball modes we refer to |lCi|; the original work was done by Keller- Rubinow. 

A new feature in comparison to the disc is the existence of Gaussian beams (a bouncing 
ball mode) along the minor axis, which is a stable elliptic bouncing ball orbit. Such bouncing 
ball modes do not exist in the disc and occur when m is fixed and n ^ oo. 

The eigenfunctions which are the Gaussian beams are characterized as follows: Since the 
minor axis is 

I = {ip,ip) e [0, pmax] X [0, 2n]; ip = -}. 

one looks for eigenfunctions with mass concentrated along this interval. Consider the ex- 
tremal energy levels that satisfy 

c^X-^ = 1 + 0{X-^). 

One rewrites (12.11) in the form: 

-X'^d^G + {X~^c^ cos^ ip-l)G = 

(2.2) 

-X-^d^F + {X^^c^ cosh^ p + 1)F = 

Given the choice of energy level, 

X-^c^ cos^ <f-l = cos^ + 0{X~^) = - sin2(y5) + 0[X'^). 

The potential ^(v?) = cos^ Lp — \ has a nondegenerate minimum at (y9 = |. So, the solutions 
G = Gm,n to the first equation in (12. 2p are asymptotic to ground state Hermite functions. 
More precisely, 

GmAv; X) = c^,„(A)e-^'=°^''^(l + 0{X-')). (2.3) 

In the second equation in (12. 2p the potential is V2(p) = cosh^ p + 1 + 0{X^^) > for A > Aq 
sufficiently large. In this case the solution has purely oscillatory asymptotics: 

^m,n(p; A) = e^^/oVcosh2x+ldx^_^^^. ^ ^2.4) 

where a-|-(p;A) ~ S^lo '^±j(/')A determined by the Dirichlet or Neumann boundary 

conditions. Moreover, from the L^-normalization condition Jj |\E'm^„(p, f )pc?p = 1 it follows 
that Cm_„(A) ~ X^^^. 

From (12.31) and (12. 4p . the Gaussian beams are roughly asymptotic to superpositions of 
g±4fcsg-Ai/ |BBj ). where s denotes arc- length along the bouncing ball orbit and y denotes 
the Fermi normal coordinate. It follows that outside a tube of any given radius e > 0, the 
Gaussian beam decays on the order 0(6"^"^ ). Hence on any curve G which is disjoint from 
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the bouncing ball orbit, the restriction of the Gaussian beam to C saturates the description 
of a 'good' analytic curve. 

2.3. Remarks, (i) The goodness requirement (11. 4p on any interior curve C G Q is implied 
by an exponential growth estimate involving only the Cauchy data {'^x\c,d^f\\c) along C. 
This is a consequence of the following unique continuation argument. 

Assume that C is a closed curve in the interior of Q. Let Uc be the domain with boundary 
C U dQ. It follows from the Sobolev restriction theorem that 

Let int{C) be the interior of the domain bounded by the curve C and take x G int{C). From 
the potential layer formula (see l3.4p (px{x) = J^{d^^g)G{x,q; \)(px{q)~G{x,q; X)dygipx{q))da{q) 
and so, by squaring both sides, using the bounds \G{x,q,X)\ = 0{1) and \dy^G{x,q; \)\ = 
(9(A^/^) and applying Cauchy Schwartz, one gets 

yxWhitntiO) < CX{yx\\l2(^c) + \\duVx\\h{c))- (2-6) 
By a standard Carleman estimate/unique continuation argument [EZt ITa[ ITa2] : 

||</3a||^i/2(c;^) < e^^yx\\l2^i„t^c)) < ^e^^(llv^A||i2(c) + Wdu^xWh^c))^ 

where the last inequality follows from (12.61) . Inserting the last bound on the RHS in (12.51) 
yields the comparison estimate relating Cauchy data along C and dQ: 

W^xh^an) < {\\v\\\lhc) + \\duV\\\L^(c))- (2.7) 
As an immediate consequence of (12. 7p we note that (II. 4p follows from the exponential bound 

\\dv^x\\L-^(C) < (i^^\W\\\L^{c) (2.8) 
involving only Cauchy data along C. 

A natural question is whether (12. 8p is automatically satisfied when ipx does not vanish 
identically on C? We hope to address this in future work. 

(ii) In the case of ellipses, it is elementary to obtain the analytic continuations of the 
boundary values of Neumann eigenfunctions to the complexification of the boundary ellipse. 
One sees that the growth rate of these analytic continuations are determined by the angular 
momenta of the eigenfunctions, i.e. the eigenvalue of the boundary Laplacian, and not by 
the interior eigenvalue. The estimates of this article are however in terms of the interior 
eigenvalue. It would be interesting to strengthen the estimates of this paper to obtain 
precise estimates on the growth rate of analytically continuations of boundary values of 
eigenfunctions. 

3. HOLOMORPHIC EXTENSIONS OF EIGENFUNCTIONS TO GRAUERT TUBES 

It is classical that solutions of the Helmholtz equation on a Euclidean domain are real 
analytic in the interior and hence their restrictions to interior real analytic curves admit 
holomorphic extensions to the complexification of the curves. A classical presentation can 
be found in [G] §5.2, and to some extent we try to conform to its notation; see also |MNt |Sl 
ILew2^ [Mort W\. Moreover, the Cauchy data along the boundary of eigenfunctions satisfying 
Dirichlet or Neumann boundary conditions admit analytic continuations into a uniform tube 
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in the complexification of the boundary, independent of the eigenvalue. We now use integral 
representations for the analytic continuations to obtain upper bounds for the growth rate of 
the complexified Cauchy data of eigenfunctions in a fixed Grauert tube of interior curves or 
curves on the boundary as A — >^ oo. 

3.1. Complexification of domains Q and analytic curves C G Q. When (M, g) is a real 
analytic Riemannian manifold without boundary, then there exists a complexification Mq of 
M as an (open) complex manifold, and Laplace eigenfunctions extend to a maximal tube in 
Mc- This has been studied in |Boul IGSll IGS21 IGLSi E] using the analytic continuation of 
the geodesic flow and wave group. The notion of maximal tube depends on the construction 
of a special plurisubharmonic exhaustion function p adapted to the metric g in the sense 
that iddp restricts to the totally real submanifold M C Mq as g. The function ^/p{t) is the 
distance from t to i. 

The analogous results for manifolds with boundary have not apparently been studied 
before. In this section, we study the complexification of the boundary and analytic contin- 
uations of the Cauchy data of eigenfunctions on the boundary (or on interior curves) for 
domains f2 C with the Euclidean metric. The complexification of the ambient space is of 
course and its Grauert tube function is \'^t\. The novel features concern the influence of 
the boundary on analytic extensions of eigenfunctions. 

We adopt the following notation from Garabedian [G] and Millar |Mll IM2] : We denote 
points in and also in by (x, y). We further write z = x + iy,z* = x — iy. Note that 
z, z* are independent holomorphic coordinates on and are characteristic coordinates for 
the Laplacian, in that the Laplacian analytically extends to ^^pr- When the boundary is 
real analytic, or when we are dealing with a closed analytic curve C C fi, we parametrize C 
by a real analytic parameterization Q : — > ~ C. We complexify C by holomorphically 
extending the parametrization to Q'^ on an annulus 



for e > small enough. Note that Q extends holomorphically to A{e) as Q*. Throughout 
the paper, the subscript c or superscript denotes the holomorphic continuation of a curve 
or function; sometimes we omit the sub or superscript for notational simplicity. 
We also use the notation 



to write the parametrization as a periodic function on [0, 27r] and often put q{s) = qi{s) + 
^'?2(s), ^(s) = qi{s) — iq2{s). This parametrization analytically continues to a periodic func- 
tion on [0, 27t] + i[—e, e]. Both s and t can be either real or complex. 
We denote the boundary data of the eigenfunction by 

ux{s) = ipxiqi{s),q2is)) (Neumann); Uxis) = d,,(fxiqiis), q2is)) (Dirichlet), 

and again write -u^ for its holomorphic extension. 

Next, we put r'^ = ~ x)"^ + ij] ~ vY so that for s G M and t G C, 





q{t) = g(e^*) 



(3.2) 




(3.3) 
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We denote by -f- the not-necessarily-unit normal derivative in the direction iq'{s). Thus, in 



dn 

terms of the notation 4- above, 4- 



d , 1 
— log r = - 

ds ^ 2 



q\s) 



+ 



Q{s) - q{t) q{s) - q*{t) 



\q\s)\^. One has 



d . —i 



q'{s) 



(l{s)-q{t) q{s)-q*{t) 



When we are using an arc-length parameterization, 



_ _ d_ 

dn dv ' 

To clarify the notation, we consider the case of = dfl. Then, q{s) = e^^ , t = 9 + i^, 



q{e + tO = e 



-i{e+i0^q* = ei(e-iC)^ and 



4 sin 



2 (g-s+ig) 
2 



Thus, logr^ = log(4sin2 ^^^ft^). Clearly, ^ = ^, so 



d , 2 
— logr 

ds 



le 



+ 



-26 



d , —i 



-le 



3.2. Layer potential representations. The representations we use are analytic continua- 
tions of layer potential representations of solutions. Let G{X,x,y) be any 'Green's function' 
for the Helmholtz equation on Q, i.e. a solution of (— A — A^)G(A, x, y) = 5x{y) with x,y & Cl. 
We will always use restrictions of the global Euclidean Green's functions on M^. 
For any closed curve C G Cl bounding a domain int(C) Green's formula gives 



{duG{X, X, q)(px{q) - G{\ X, q)dyipx{q)) da{q), 



(3.4) 



c 



an dva 



where da is arc-length and where dy is the normal derivative by the interior unit normal. 
We will analytically continue this formula for various choices of analytic C. 

When C = dQ and the eigenf unctions satisfy standard boundary conditions, the formula 
simplifies. In the case of Neumann eigenfunctions (px in Q, which are emphasized here, 

d 

-G{\,x,q)ux{q)d(7{q), x G (Neumann). (3.5) 

Here, and henceforth, we denote the restriction of ip to dQ by uj. In the Dirichlet case, the 
corresponding formula is 

f ~ 9 

ipx{x) = - G{\,x,q)—ipx{q)d(T{q), (Dirichlet) (3.6) 

Jon OUq 

To obtain concrete representations we need to choose G, and we often choose the real 
ambient Euclidean Green's function S (in the notation of [GJ, §5), 

S{\x,y-i,^) = -Y^'\\T), (3.7) 

where r = \/zz* is the distance function and where 



y(o)(Ar) = Jo(Ar)log(fcr) - ^ 



(-1)^(1 + 2- 



m=l 



(3.^ 



is the Bessel function of order zero of the second kind. The Euclidean Green's function has 
the form 

1 



S{\ v; X, y) = A{X, ^, T]; X, y) log - + 5(A, r/; x, y), 

r 



(3.9) 
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with 



( 1 \k (Ar) 



A = Jo(Ar) := Er=o(-l 



2fe 



22fe(fe!)2 ■ 



^ - - 2^m=l 4'"(m!)2 ^ -'ol^r j lOg A. 

The coefficient A is known as the Riemann function and B is the Bessel function of order 
zero of the ffist kind. 

Sometimes it is more convenient to use the (complex valued) Euclidean outgoing Green's 
function ¥Lal^\kz), where Hag^'* = Jq + iY^ is the Hankel function of order zero. It has the 

same form (13.91) and only differs by the addition of the even entire function Jq to the B term. 

a 



By the 'jumps' formulae, the double layer potential ■SzS{X, x, q) restricts to the boundary 



as \5q{q) + -^_S{\,q,q) (see e.g. [T]). Hence in the Neumann case the boundary values 
satisfy, 

[ d 

ux{q) = 2 —S{\,q,q)ux{q)da{q) (Neumann). (3.10) 
Jan 

In the Dirichlet case, one takes the normal derivative of ipj at the boundary to get a similar 
formula for dy(pj\dQ,, with a sign change on the right side. We have, 

^5(A, g, ~q) = -XY}'^ (Ar) cos Z(g - q, z/,-) (3.11) 

where 

nY.iz) = ^ + 2Mz) (log(^/2) + 7) - EZii-^^' InizTy. (^^''-'[1 + 2 ELi 

(3.12) 

Here, 7 is Euler's constant. As is well-known, the pole of Y^ is cancelled by the cos Z{q—q, Ug) 
factor. 

If instead we use the Hankel free outgoing Green's function, then in place of (13. lip we get 
the kernel 

NiX,qis),q{s')) = Ha(^)(A|g(s) - y|)|,=,(,,) 

(3.13) 

= -AHai^)(A|g(s)-g(s')|)cosZ(g(s')-g(s),j/,(,,)). 

We now consider the analytic continuations of these formulae. First we parametrize C 
by a real analytic parameterization. Since Jq is even, A{C,,T],x,y) admits the holomorphic 
continuation R{C, (*, z, z*). To simplify notation put R{q{s), q{s), z, z*) := R{s; z, z*). 

3.3. Interior curves. First we consider the simple problem of analytically continuing the 
representations (13.51) and (13. 6p . We are interested in restrictions to real analytic closed curves 
C C 

Let q : [0, 2tt] —* C denote a real-analytic parametrization of the interior curve C. For 
e > sufficiently small, we consider the annulus A{e) = {t G C; < \t\ < e^} and the 
corresponding complexification of C given by 

Cc = {Q^(e*(^+^")) G g^(A(e)); Q^{e'^'+'^^) = q^{s + tr); < s < 27r, -e < r < e}, 

where, we recall that denotes holomorphic continuation. Since C is assumed to be an in- 
terior curve, it follows by compactness of dQ that for \'^q'^\ sufficiently small, r'^\ccxan 7^ 0. 
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As a result, one can choose a globally defined holomorphic branch for logr and so, the holo- 
morphic continuation formula for Neumann eigenf unctions in this case follows immediately 
from (13. 5p : 

^l{q''{t))= ! ^SiX,q''{t),q{s))u,{s)da{s). (3.14) 
Jan duq 

3.4. The case C = dQ. We let q{s) denote a real analytic paramaterization of dQ. We use 
the arc-length parameterization so that da{s) = ds. From (13.91) . we can write (13.101) as 

«A(t) = j-J-',{-u,{s)^^{s,t))\ogr^ds 

(3.15) 

-i />,(.) A(., t)ifrf. - i /;(-«,(.)f(., t)ci.. 

With different choices of B the same formula is valid for the outgoing Green's function as 
well. 

Since the integral is now over dVt, the logarithmic factor in S gives rise to a multi-valued 
integrand, and it is not obvious that the representation can be holomorphically extended. 
We first concentrate on the case where C = dfl when dQ is real analytic, and later consider 
the case where C is an analytic arc of dQ when it is piecewise real analytic. 

So assume at first that C = dfl is real analytic, and as above let u\{s) = (px{q{s)) be the 
boundary traces of the Neumann eigenfunctions. Our goal is to analytically continue the 
representation (I3.15p . 

Millar's formula for the holomorphic continuation of the Cauchy data is as follows; let 

^t;z,z*)= / ux{s)—R{s,z,z*)ds. (3.16) 

Proposition 7. |Mll IM2j The boundary data u = ux of the eigenfunctions of the Neu- 
mann problem admit the following holomorphic extension to a uniform tube around dQ in 
its complexification {dQ)c: (for '^t > 0, < 0) 

u^.it) = ±t^t, g(t), q*it)) + Ms; g(t), g*(t)) + ^Uxis)Ris, g(t), g*(t))] ^£^ds 
+ Ms; q{t), m) - ^Ms)R{s, g(t), q*m ^^f^ds 
-2j^ux{s)^{s;q{t),q*{t))ds. 

(3.17) 

Proof. Since we depend crucially on this Proposition, and since it does not appear to be 
well-known, or to be proved in detail in |Mll IM21 IV] . we supply the details of the proof. 

We will analytically continue the formula (13.151) . Although u is real analytic on dQ and 
hence admits an analytic continuation to a small complex 'tube' {dQ)c, it is not clear that 
the representation (I3.15P can be extended analytically due to singularities of the integrand. 
Moreover it is not clear that the right side of (13.270 is in fact complex analytic. The main 
task in the proof is to clarify these points. 

We begin by showing that the last two terms of (I3.15P analytically continue in a straight- 
forward way. 
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Lemma 8. The integrals (i) ^ uxis)A{s,t)l^{s,t)ds, resp. (ii) ^ ux{s)^{s,t)ds, 
are real analytic on the parameter interval parametrizing dD, and are holomorphically 
extended to an annulus by the formulae 

d 



resp. 



u,{s)^{s;q(t),q*(t))ds. 

Proof. Any derivative of log is unambiguously defined and we have 

1 dr dlogr 1. q'{s) q'{s) , 

r dn dn 2i q{s) — q{t) q{s)—q*{t) 

In the real domain, q*{t) — q{t), so 

1 dr (/{s) 



r dn q{s) — q{t) 
We recall that |^ = \q'{s)\~^^. In terms of the real parametrization q{(p), 

dr 

— = cos Z{q{if 2) - 

vanishes to order one on the diagonal in the real domain so that is real and continuous. 
In complex notation, the same statement follows from the fact that 

t-^s s-t q{s)-q{t) s-t ^ 

where ^ is real when s,t G M. Hence, ^^(^^^(Tj is continuous for s, t G [0,£] and since 

q{s),q{t) are real analytic, the map s — >• [ g(g)ig^(t) ~ q{s)-q*{t) ^ ^ continuous map from 
s G [0, £] to the space of holomorphic functions of t. 

Since A — Jo(kr) is an analytic function of r^, A^^^ has the form 

1 r I'i^) I'i^) 



2i'q{s)-q{t) q{s)-q*{t)'' 

for an analytic function F. Clearly, F(r^(s, t)) is also a continuous map from s G [0, i] to the 
space of holomorphic functions of t. Hence, so is the product and therefore so is the integral 
over s G [0, i] of the product. 

Similarly for case (ii). In this case, B is an entire function H{r^) of which is of the form 
r'^h{r'^) for another entire h. Hence, 

- r^H'ir^) X ^ ^ ^'^'^ ^ 



du 2i q{s) — q{t) q{s) — q*{t) 

So the integral (ii) also admits an analytic continuation. 



□ 
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Thus, the difficulty in analytic continuation of the representation is entirely with the 
integral J^{— u{s)^{s,t)) log r'^ds. Due to the logarithm, the analytic continuation of the 
integrand is multi-valued in any neighborhood of dQ. Nevertheless, the integral admits a 
single-valued analytic continuation there. 

Lemma 9. The integral J^{— u{s)^{s,t)) log r'^ds extends to a holomorphic function oft in 
a neighborhood of dQ in {dQ)c given by 

«(t, q*{t)) + 1^ [$(.; ^^^^^ds + j' [$(.; g*(t))] ^(^p^^^, 

where, ± corresponds to Qt > 0, < 0. 



Proof. We first observe that 



W7r = Mr)^- (3.18) 



du du dp 

Now J I is odd in r so fl3.18p has the form 

_ 77/ 2 "^logr _ 2x 2 1 r Q-'i^) q' { 



Fir') r- ^ F(r') H ^ ^ F(r') r' [^l^^^ " ^^ff^l j . (3.19) 

where F is a holomorphic function. From fl3.19p it follows that ^{t, s) is a smoothly varying 
family of holomorphic functions of t in a sufficiently small annulus. 

Thus, our problem is a special case of the general problem of analytically continuing the 
integral f{s) log r'^{s,t)ds where / is real analytic and where r'^{s,t) is given by (13. 3p . 
In our case, f{s) also depends holomorphically on t but this does not affect the analytic 
continuation issue. 

To define the analytic continuation, we specify a branch L{s, t) of the multi-valued analytic 
continuation of logr^(s, t). We ffist slit the complex parameter annulus through the vertical 
segment through to obtain the complex t parameter strip / = [0, i] + i{—e, e). Our integrals 
only involve pairs {s, t) G [0, i] x /. We then remove the the set < s < ^t. For fixed s, these 
cuts disconnect the t-annulus into two strips. On one strip, < 3ftt < s in the parameter 
interval, while s < < £ in the other. We then further slit these strips along the real 
segment 9t = of /, to obtain four strips or 'quadrants'. In the right 'half-plane' where 
s > JJt, we define '^L{s,t) so that it is continuous in the right 'half-plane' and tends as 
Sjt — s> from either top or bottom. In the slit left 'half-plane', {s < 3?t}\[0,3?t] we define 
L{s,t) by continuation from the right half plane. It then tends to ^2tt as 53t ^ from 
above, resp. below the cut along the real axis s < 3?t. 

To illustrate, we consider the basic case of the circle, where q{t) = e** and where we 
are defining arg((e** — e**)(e~*'' — e~**)). We fix = to and consider the map (s, r) 
^^is _ g«i^(g-«s _ g-it-) -^here t = to + it. In the 'ffist quadrant' s > to, '^t > 0, this map is 
anti-holmorophic and takes the real axis Q't = to the positive real axis and the 'imaginary 
axis', s = to and Q't > 0, to the negative real axis. Since the map is anti-holomorphic, 
the image of a counter-clockwise path in the first quadrant from the real to imaginary axis 
is a clock-wise path from the positive real axis to the negative real axis, so the arg equals 
—71 on the imaginary axis. As the path in the domain moves counter-clockwise in the 
second quadrant to s < to, ^t = the image path moves to argument — 27r. Similarly, the 
continuation in the fourth and third quadrants leads to a value of 27i on the axis s < to- 
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The following claim is the key one to obtain a single valued analytic continuation (cf. 
Millar [M2] l 

Claim: If / admits an analytic continuation to a neighborhood of dfl, then the integral 
Jq fis) \ogr'^{s,t)ds admits an analytic continuation to a neighborhood of dil in (9f2)c by 

/ f{s)\ogr^{s,t)ds / f{s)L{s,t)ds±2TTi / f{s)ds, (3.20) 
Jo Jo Jo 

where the path from to t is defined in the integral is the same as the path used to analytically 

continue logr^(s,t), and where the + sign is taken for > and the — sign when < 0. 

First, we check the periodicity of the right side in (13.201) . This follows from the fact that 

pi nl+i^t pi 

/ f{s)L{s,£ + m)ds±2TTi / fis)ds= / L{s,m)f{s)ds 
Jo Jo Jo 

ri+i'^t 



±2m / f{s)dsT'2m / f{s)ds 
Jo Jo 



L{s,m)f{s)ds±27ri / f{s)ds, 

Jo 

where, the last line follows from the Cauchy formula and the fact that f{t + i) = f{t). 
To prove the claim it suffices to show that the right side in (13.201) is 

• (i) Holomorphic in the upper annulus Q't > 0; 

• (ii) Holomorphic in the lower annulus St < 0; 

• (iii) Continuous in the whole annulus, and restricts to Jg f{s) \ogr'^{s,t)ds for real t. 
Let us first recall why this is sufficient: 

Lemma 10. (see e.g. |Mu] ) Let Q (Z C be a domain, and suppose that C dVt is a closed curve 
separating VL into two domains D~^, D~ with common boundary C. Suppose that F~^, F~ are 
holomorphic on , resp. D~ and that = F" on C . Then the function F defined by 
F = F+\ B+-,F I/)- and F^ on C is holomoprhic on all ofQ. In other words, a sectionally 
holomorphic function which is continuous is holomorphic. 

Let us prove (iii) first, since it explains the second term on the right side of Claim 13.41 
With no loss of generality, suppose that 53t — O"*" with t ^ t^. Then 

pt p£ pto 

f{s)L{s,t)ds + 2m / f{s)ds^ / f{s)L{s,to)ds + 2m / f{s)ds, 



JO Jo 
and we must show that 

„2/ 



f{s)L{s,to)ds + 2m / f{s)ds= / f (s) log r {s, to) ds. 
Jo Jo 

Here, argr^(t, s) = while '^L{s,t) equals zero for s > t and equals —2it for s <t. Hence, 
the imaginary part of the left side cancels and we obtain the right side. 

Now let us prove (i)-(ii). Since the proofs are essentially the same we only prove (i). 

We first note that all branches of analytic continuation of log r^(s, t) differ by constants in 
2'KiZ. Hence, if the period (/) := j f{s)ds of / vanishes, then all choices of branch of log "^^ 



r 



give the same value of the integral f{s) log r'^{s,t)ds. Similarly, the integral f{s)ds is 
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only multi-valued due to the period of /. Hence, when the (/) = 0, both terms on the right 
side of the claim are well-defined independently of any choice of integration path or branch 
of logr^. Since we can write / = (/ — (/)) + (/), we only need to show: 

(1) f{s)L{s, t)ds ± 2Tii f{s)ds is holomorphic for > when (/) = 0; 

(2) £ L{s,t)ds±2m j^ds is single- valued holmorphic function for 53t > 0. 

To prove (1), we assume (/) = and let F{t) = f{s)ds be the (well-defined) primitive 
of F in the annulus. We then integrate by parts in the first integral to obtain 

F'{s)L{s,t)ds 

= Fii)Lii, t) - F(0)L(0, t) - Fis)^^ds - Fis)^^ds (3.21) 

= -i:Fis)^ds - J^Fis)^^ds 

Here, we use that F{i) = F{0) since it is a single-valued holomorphic function on the annulus 
and that F{0) = by definition. It is clear that the right side of (13.211) is holomorphic in 
53t > since poles occur only when = 0. Since F{t) is single valued and holomorphic, 
this proves (1). 

To prove (2) we write 

j:Lis,t)ds±2mJ^ds = C\og ^'^-'^r''''f2-'^T' ds 

(3.22) 

+ J^^ log ((e^^ - e^*)(e-^" - e"^*)) ds ± 27ri ds, 

We observe that the first term is holomorphic for > since the arg of both numerator and 
denominator are continued so that each arg tends to —27r as St — for s G [0, 3?t] and so 
that each arg tends to zero for s G [3?t, i]. Hence, the arg of the ratio tends to zero as 53t — 
in both integrals. Since the arg is only ambiguous up to a constant in 27rzZ, it follows that 
the arg of the ratio is well defined and single valued and the integrand is well-defined as a 
single- valued holomorphic function for > 0. Therefore, to prove (2) it suffices to show 
that for '^t>0 

g{t) := / L{s,t)ds + 2ni ds = 0, (3.23) 
Jo Jo 

where L{s, t) = log ((e*'^ — e**)(e~*'' — e~**)) with our choice above of the branch cut at s = 3?t. 
Note that g is an analytic continuation of the integral J^^ log |e** — e**p(is = for real t, so 
analyticity of g is equivalent to g = 0. 
This reduces the analysis to the integral 

/'27T /'27r 

/ log((e*"-e**)(e'^^-e-**))cis= / log (2 - 2 cos(s - t)) ds, 
Jo Jo 

where as above the logarithm is defined by breaking up the integral into jj** + J^^ and 
defining the arg by the method above. Note that formally the integral is constant in t by a 
change of variables but that such a change of variables is not consistent with the definition of 
the logarithm. However, the integrand is a function of s — t and so, 4 log (2 — 2 cos(s — t)) = 
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-;jjjZs) log (2 — 2 cos(t — s)) is well-defined independent of the branch of log used (since these 
differ by integer multiples of 27ri). Hence, 

|/o"log(2-2cos(s-t))ds = log(2-2cos(s-t))ds 

= -log(2 - 2cos(s -t)) 



|27r 

lo 



= —2711, 

by definition of the logarithm. It follows from (13.231) that 
d d f'^'" d 

—g{t) = — log (2 - 2 cos(s - t)) ds + —{2mt) = -2ni + 2m = 0. 

CLT db J Q db 

Hence g is constant and as noted above it equals for real t. 

This completes the proof of the Claim and hence of the Proposition. 



□ 



□ 



Remark: By integrating by parts directly in the integral Lf{t) := logr^(s, t)/(s)(is for t 



real and using that J^^log \e^^ — e^^\'^ds = 0, one gets the formula 

(3.24) 

where, F{t) := /q(/ - {f))ds. It follows from (13:241) that for t G [0,i] + i[-e,e] there is an 
alternative formula for the analytic continuation of Lf which is given by 

^'^^'^^^ - - II ^ ^(Iftm) ■ - ''''' 



We note that the log in the second term on the RHS of (13.251) is defined unambiguously 
(independent of branch) since as St ^ and s — >■ from either side, we have that 
aig[Q{e^^) — (5(e**)] — arg[e** — e**] and so the arguments cancel. The same thing is 
true for the ratio involving Q*. 

Remark: It should be noted that the proof only makes use of the fact that ux admits a 
single- valued analytic continuation to {dfl)c- The right side of (I3.15P defines an operator Nu 
on C{dfl), and the proof shows that although Nc{X, (, q) is multi- valued, its integral against 
u admits a single-valued holomorphic extension as long as u admits one. Let us check the 
argument in the case of the unit disc, where is a convolution operator with kernel 

iV(A,^,^) = ^AHa«(A)j;(A)e-(^-^). 

nez 

The Fourier coefficients A Ha^^-'(A) J^(A) have only a slow decrease reflecting the singularity 
of log I sin^(^ ~ V^)|) so the kernel does not admit a holomorphic extension e*"^ — > in a 
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pointwise sense. However, its integral Nu against a real analytic function u ~ Y2n '^nC*"''' ad- 
mits the holomorphic extension J2n&i. ^^^n\^) '^nW^nC"' to the maximal domain to which 
u itself holomorphically extends. 



3.5. The case C fl dQ ^ with C ^ dQ. In this case, we cannot use the boundary 
conditions to simplify the integral 

ifxiq) = 2 (^J^^_S{q,q;X)ipx{q)da{q) - J^S{q,q;X)^_ipx{q)da{q)^ . (3.26) 

The first term is handled exactly as in the case C = dfl, while the second term (the single 
layer potential term) is new. 

To continue the second term we write it in the form 

-i / A{t,s;X) \ogr'^{s,t)du'Px{s)da{s) - [ B{t, s; X) d^^^xis) da{s). 
^ Jc Jc 

Since B = F{r^) where F is entire, the second term above has a straightforward analytic 
continuation. The first term is another case of the integrals discussed in the previous section, 
and its holomorphic continuation is: 

R{q{t), q*{t)- s)L{s, t)d,^x{s) da{s) T 27rz / R{q{t), q\t)- s)d,{s) da{s). 

Jo 




3.6. The case where C is piecewise real analytic. We now consider the case where C 
is piecewise real-analytic. This has previously been discussed in |M3] . 

By a piecewise analytic embedded closed curve C of length L{C) we mean a curve of the 
form C = IJj=i where 

• Cj C are the maximal real analytic components of C, enumerated in counter- 
clockwise order so that Cj intersects only Cj-i and Cj+i. 

• The Cj are parameterized by m real analytic functions qj(tj) : [0,ij] Cj on m 
parametrizing intervals (where £j = L{Cj) is the length of Cj. We assume CjflCj-i = 
{qA^) = Ij-ii.^j)} when m > 2. 

We denote the Cauchy data of the eigenfunction ipx on the boundary component Cj by 
u\. Our aim is to analytically continue u\ to Uj=i ^i,c where Cj^c is a complexification of 
the interior of Cj. Thus, as we define it, Cc is pinched at the corner points and the analytic 
continuation of the boundary data of (fx is somewhat simpler than in the fully analytic case 
in that we are analytically continuing to a smaller kind of neighborhood of C. 

Millar's formula for the analytic extension of u\ to in the Neumann case is given by: 

Proposition 11. [M3] The boundary datau\ of the eigenf unctions of the Neumann problem 
admit the following holomorphic extension to a uniform tube around the interior Cj of Cj 
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in its complexification: (for '^tj > 0,<0) 

(3.27) 

+ /o M^n; qit,), qit,)) - mA(s„)i?(s„, ^(t,))] ^^£z^)ds 

-2 /o MA(s)|f q{tj), q*{tj))dsn. 
Proof. We only sketch the proof, because it only involves a small modification of Proposition 

El 

First, Green's formulae fl3.4p - fl3.5l) remain correct in the piecewise analytic case, with the 
definition that on the jth component, the normal derivative is calculated by taking the limit 
from within the jth component. 

The verification of the Millar formula is then similar to the fully analytic case. The 
main difference is that we now have pairs {sn,tj) of parameter points which may come from 
different intervals. When n = j there is no difference in the argument except that Cjc is 
not an annulus but rather two regions meeting along a common interval. But the same 
choice of branch of the logarithm extends Uj holomorphically above and below the interval, 
and the first term on the right side ensures that the two holomorphic extensions agree on 
the common interval. When n ^ j, one defines argr^(s„,tj) = for all real Sn,tj. Since 
Qnisn) 7^ Ijiij) when n ^ j the logarithm extends to a holomorphic function in tj with this 
choice of branch. 

□ 



4. Growth of zeros and growth of 

Let Cc be the complexification of a real analytic curve C dVt. The purpose of this section 
is to give an upper bound for the number of zeros of in the annulur region q^{A{e)) where 
A{e) = {z G C; e"'' < \z\ < e^}. For \j G ^^(v^) and for a region D C Cq denote by 



n{X,,D) = #{g^(t) G D : ^^g^W) = 0}. (4.1) 

The following estimate is suggested by Lemma 6.1 of Donnelly-Fefferman |DF] . 

Proposition 12. Suppose that C is a good real analytic curve in the sense of (T^. Nor- 
malize U\ so that ||ma||l2(c) = 1- Then, there exists a constant C{e) > such that for any 
e > 0, 

n{X,Q^{A{e/2)))<C{e) max log|«^(g^(t))|. 

g^{t)GQC(A(£)) 



Proof. Let denote the Dirichlet Green's function of the 'annulus' Q'^{A{e)). Also, let 
{a^ jjS'^''^^^'/^^" denote the zeros of in the sub-annulus Q^{A{e/2)). Let Ux = p^jp^ 
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where | |qc(a(,)) = m&x^^QC(^A(e)) \u{C)\ - Then, 

log\U,{q^m = /Qe((^(,/2)))Ge(g^(t),«;)aaiog|«,(«;)| + i/,(gC(t)) 

since dd\og\u\{w)\ — '^a eCc-u'^{a )=o^aj- Moreover, the function Hx is sub-harmonic on 
Q^iA{e)) since 

ddHx = ddlog\U,{q''m- E ddG,{q''{t),a,)= J] > 0. 

ajGQC{A{e/2)):ux(a,)=0 ajeQ'^{A{e))\QC{A{e/2)) 

So, by the maximum principle for subharmonic functions, 

max Hxiq'^it)) < max Hx(q^(t)) = max log\Ux(q^(t))\ = 0. 

QHM^)) 9QC(A(6)) 

It follows that 

log|C/A(g^W)|< E G,(g^(t),a,), (4.2) 

ajeQHA{e/2)):ux{aj)=0 

hence that 

T^.^ loglUxiq'^m < ( max G,{z,w)) n{\ {A{e / 2))) . (4.3) 

Now Ge{z,w) < max^gQC(a^(^)) Ge(2;, -u;) = and Ge{z,w) < for w G Q'^{A{e/2)). It 
follows that there exists a constant z/(e) < so that niax2^gQC(^(^/2)) G^e(-S;'ii') ^ ^(^)- Hence, 

max log|C/A(g^(i))| < i.(6) n(A,g^(^(6/2))). (4.4) 

Since both sides are negative, we obtain 

n(A,QC(^(e/2))) < ^|niax,c(,)eQC(^(,/2))log|t/A(g^(t))|| 

^ ^ (maXqC(t)6QC(^(,)) log |xiA(g^(i))| - maXgC(t)gQC(^(,/2)) log \ux{q^{t))\) 
^ WT\ maXgC(t)gQC(^(,))log|ifA(g'^(^))|, 

(4.5) 

where in the last step we use that inaXqC(t)gQC(^(£/2)) log |wA(9'^(i))| > 0, which holds since 
\ux\ > 1 at some point in Q'^{A{e/2)). Indeed, by our normalization, ||wa||l2(c) = 1; and 
so there must already exist points on the real curve C with \ux\ > 1- Putting C{e) — j^^^ 
finishes the proof. □ 

Remark: An alternative approach is to use Jensen's formula, 

/XA, qc{A{p))dp = Mf,{e) + Mf,{-e) - 2My,(0), (4.6) 

where 

^ ^ \de 



MAp)^(^[ log|/A| 

\27T 
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However, this would require an analysis of the real logarithmic integral Mj^(O). As the 
example of the Gaussian beam shows, Mj^(O) may be of order A due to exponential decay 
away from the 'classically allowed region'. We plan to analyze such integrals elsewhere. 



5. Proof of Theorem [5] 

Proof. We prove Theorem [5] before Theorem [2] since it is easier. The proof of Theorem [5] 
uses the analytically continued potential layer formula fl3.14p to bound maxQC(^(^)) log 
from above. Then, an application of Proposition [T2] gives the result. 

5.1. Upper bounds for the analytically continued eigenfunctions. Let q : [0, 27r] C 
be an arc-length parametrization. We denote the parametrization of the boundary, 5^2, by 
q : [0, 27r] dVL. In this case, the formula for holomorphic continuation of eigenfunctions is 
given by (13.141) : 

^l{q^{t))=j N{\,q^it),qis))u,is)dais), (5.1) 
Jan 

From the basic Hankel function formula (IS.lSp for N{X, q, q) and the standard integral 
formula 

/ 2 \ ^ pi(2-37r/4) /-oo 

<'W=y -TlmL (5,2) 

one easily gets an asymptotic expansion in A of the form: 

A; 

N{q^{t), qis); A) = e^^^^'^'^W'^'^^)) J] a,{q^{t), q{s)) A^"^' + C)(e^Ar(,C(i),,-(,)) ^i-fc-i^^ ^5 3) 

j=0 

Note that the expansion in (15. 3p is valid since for interior curves, 

Co:= min |g(t) - g(s)P > 0. 

iq{t),q{s))£Cxan 

Then, ^r'^{q'~"\t),q{s)) > as long as 

l^q^'m^KCo. (5.4) 

So, the principal square root of has a well-defined holomorphic extension to the tube (15.41) 
containing C. We have denoted this square root by r in (15.31) . 

Substituting (15.30 in the analytically continued single layer potential integral formula (15. ip 
proves that for t E A{e) and A > sufficiently large, 

^^(g^(t)) = (27rA^) / e^^^(«"W'^>))ao(g^(t), g»)(l + O(A-i) ) «,(.)da(s). (5.5) 
Jan 

Taking absolute values of the integral on the RHS in (15. 5p and applying Cauchy- Schwartz 
proves 

Lemma 13. For t E [0, /] + e] and A > sufficiently large 

|</^^(g''(t))| < CiA^/^exp A ( max g?^r(g^(t), g(s))) ■ WuxUnan)- 

\q{s)&an I 
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From the pointwise upper bounds in Lemma [131 it is immediate that 

log max |<^a(9'^(^))I < C'maxA + C2logA + log||MA||L2(9C), (5-6) 

where, 

Cmax = max ^ir{q^{t),q{s)). 

(gC(t),g(s))gQC(^(,))xan 

Finally, we use that log ||MA||L2(an) = 0{\) by the assumption that C is a good curve and 
apply Proposition [T2] to get that n{\,C) = 0{X). 

□ 

6. Proof of Theorem El Zeros 

The proof of Theorem [5] is more complicated than that for interior curves because we need 
to invert the Volterra operator of Proposition [71 We recall that the analytic continuation of 
u is the solution of a Volterra equation, 

{I + K^)ux{t) = Uxit), (6.1) 

where Ux{t) has an explicit analytic continuation, where 

/"* (9/? 

Kxux{t) = —(X,s,q{t),q*{t))u,{s)ds (6.2) 
in Millar's notation. Here, R = A is the Riemann function and so explcitly, 

Kxuxit) = t ^^^{X,s,q{t)q*{t))ux{s)ds. 



du 



Therefore, 



K,{t,s) = l[o,](s)Ji(Ar)r^^(t,s) = l[o,](s)r Ji(Ar) ^ ^^'^ ^'^'^ 



du ' \<l{s) — q{t) — q{t) 

where l[o,t] is the indicator function of the interval [0, t]. We note that the pole of g(sRt_(_'^i)Lg(t)» 
at the upper limit of integration s = t is cancelled because rJi(r) begins with r^. So the 
integrand is regular and holomorphic along the path of integration. 
On the right side of the Volterra equation, 

w^(t) T ^^t, q{t), q*{t)) = [^s; q{t), q{t)) + ^Uxis)Ris, g(t), g*(t))] J^ds 



+ Jo Ms; q{t), q*{t)) - ^u,{s)R{s, q{t), g*(t))] ^£$r^ds (6.3) 



-2j^u,{s)^{s;q{t),q*{t))ds 

the Cauchy data u\ is only integrated over the real domain where by a standard Sobolev 
estimate it has polynomial growth in A. And further, the Riemann function and other special 
functions occurring there have exponential growth at most given by the ambient complexified 
distance function. The main problem is thus to invert the Volterra operator / + Kx and to 
obtain a similar growth estimate for (/ + Kx)~^{RHS). 
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We first simplify tlie operator, K\. Given t = + i'^t we may cfioose tlie contour to 
go along the real interval [0, 3?t] and then to go along the vertical line segment + is for 
s G [0, $5t]. This decomposes K\ = K^'' + where 

K^x^Mt)= / Ux{s)—RiX;s;qit),q*it))ds (6.4) 



and where 



KfKx{t) = / uxi^t + is)—R{X; m + is; q{t), q*{t))ds. 
Jo 

We move the K^^ term again to the right side since it only involves the Cauchy data on the 
real domain. 

We now write t = + i'^t and treat as a parameter. We need to study the mapping 
properties of K^^^ and (/ + Kf^)~^ on the weighted Hilbert space L^([— e, e], e~^^'^^^d^t). 

6.0. 1. Model example. As a model example, we consider the operator K^wd/) = e^^'^^^'>u{s)ds. 
A simple calculation shows that for n > 0, 

irr^(y,.) = e^(^-^)^^^, 

77/. 

and 

(/-ir,)-i(2/,3) = e(^+^)(^-^). 

Hence, in the model example, the exponential growth of the kernel (J — K\)~^{s, SJt) is the 
same as for K\{s, ^t). 

6.1. Upper bounds. In view of the growth estimate for complex zeros in Proposition [T2| 
one needs to determine asymptotic pointwise upper bounds for the |u|^(t)| as A ^ oo. In 
this section, we prove: 

Lemma 14. Given t G [0,/] e, e] and A > sufficiently large, there exists a constant 
C > such that 

\ux{t)\ < expCA|$5t| • \\ux\\L^{dn)- 

Proof. Let Co > be a constant. To bound the kernel K\ (Q't, s) we split the analysis into 
two cases: (i) \r{lR,t + is,t)\ < ^ and (ii) |r(3?t + is, t) | > 

6.1.1. The range \r\ > In this range, Ji has an asymptotic expansion given by 

k 

J^(Xr) = ^ A-3-Jaj.(r)e^^" + C(A-^-fc-ie^^"). 



i=o 



From the identity 



= l[o,^,](s)Ji(Ar)r^ 

= l[o,m] (s) r Ji (Ar) [^^^^ - ^^^^] , 

it follows that there exists a symbol S\ of order — ^ such that 

\Kf\Qt, s)\ < Sxi^t, s) l^oM^) e^l^K«*+-'*)l. (6.5) 
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The estimate (16.51) is locally uniform in and the dependence on the parameter is 
implicit . 

6.1.2. The range \r\ < In this range, the asymptotic expansion breaks down when |r| ^ 
J and so, we use the standard integral representation for Ji to get the necessary bounds for 

Since 

/•TT 

J,{Xr) = -m e'^'^"'^ cosOde (6.6) 

and I cos^l < 1, it follows immediately from (16.61) that in this range, 

|Ji(Ar)| < Ce^l^'-I, (6.7) 

and so, when \r{^t + is,t)\ < 

\Kf\^t,s)\ < Cl[o,^>*](s)e^l^^(«*+^^'*)L (6.8) 
Combining the estimates ( 16. Sp and ( 16. Sp . it follows that 

\K^x\'^t, s)\ < Cl^oM') e^l^KJJt+-,*)l (6.9) 

locally uniformly in |s| + \'^t\ and in Again, the dependence of Kx on the parameter 
has been suppressed. 

6.1.3. Pointwise estimates for r . By definition, 

r(3fJt + is, t) = {q{m + i^t) - q{m + is), q{m + i^t) - q{m + is))^ (6.10) 
Taylor expansion around s = S^t in (16.101) gives 

\r{m + is,t)\ <C\Qt- s\, (6.11) 

since, 

q{t)-q{m + is) = I ^q{m + i{tQt+{l-t)s))dt = (Qt-s) [ q'{m + i{tQt+ {1 -t)s))dt. 
Jo dt Jq 

From (16. lip and the bound (16. 9p it follows that there are constants Cj > 0;j = 1,2, such 
that 

($5t,s)| < Cil[o,st](s)e^^l^*-^l. (6.12) 
Next, we expand (/ — K^^^)^'^ in a norm convergent geometric series, 

V[irfr(S5t,3), ;[irfr(53t,s):= / / ■■■/ K^^\^t,s^)---K^^\s,,s)ds,---ds^. 

n=0 -^0 Jo Jo 

We recall that the n-simplex A„ is defined by 

{(si, . . . , s„) : < Si < S2 < • ■ ■ < s„ < 1}. 
Let ■ An be the scaled simplex. Applying the estimate (I6.12p to each factor in the above 

(2) 

formula for [K\ ]" gives the following pointwise bound: 

I [irf ^]"(53t, s)\< [ e^^(^*-^") ■ e^^(^"-^"-i) ■ ■ ■ e^^(^2-«i) . ^cxis^-s) . . . ^^^^ 

J^t-A„ 
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So, by the model example, 

To complete the proof of proof of Lemma [TU we note that for q'^{t) G Q'^{A{e)), the right side 
of the analytic continuation formula (16.31) together with the K'^^ term satisfies the estimate 

(**) < Ciexp (A max ^ir{t,s)] ■ \\ux\\mdn) < Cie^^^l^*!, (6.13) 

since by our normalization, ||MA||L2(gn) = 1. It follows that 

{I - Kf^)-\**) < C / e(^+^)^(^*-^)e^2^^(is <Cexp(Amax{C + l,C2}|5>t|). (6.14) 

Jo 

This finishes the proof of Lemma [TH □ 

So, from Lemma [m logmaXgC(()gQC(^(^)) I'/'aI'^)! — ^^^'^ Proposition [12] then implies 
that n(A;(91]) = C»(A). □ 



7. Critical POINTS: Proof of Theorems [3] and [5] (2) 

We now prove part (2) of Theorem [5] concerning the growth of critical points. It immedi- 
ately implies Theorem [3l The argument is similar to that for counting zeros, the only change 
being that we now take the derivative and restrict to the boundary in (13.61) . For the sake of 
brevity we only sketch the proof. 

In the Dirichlet case, the jumps formula for the double layer potential gives (I3.15P ex- 
cept that now u\. denotes the restriction to the boundary of We refer to [H IHZj for 
background. We then define n{Xj,D) as in (14.11) but for the new u\.. The layer potential 
representation implies the analogue of Lemma [T3] and by Proposition [T^] we conclude that 
the number of complex zeros (hence real zeros) is 0(A). 

In the Neumann case, we must take the tangential derivative u'y^{t). Since the normal 
derivative is zero, the critical points of the tangential derivative are critical points of the 
eigenfunction along the boundary. The tangential derivative now has the representation, 

+ /o(-^A(.)f (.,t))|logr2rf. (7.1) 

-^/o«A(3)(|(A(.,t)i|5))rf.-i/;(-...(.)gS(.,t)rf.. 

The analytic continuation of Proposition [7] applies equally to the equation (17. ip and the 
analytic continuation of u'{t) has the form (16. ip with a slight change in K\. However, the 
phase function is the same, so the proof of Lemma \\M applies with only small modifications 
to the new Volterra operator, giving the upper bound 

\dtu'i{t)\ < AexpCA|$5t| ■ \\ux\\L\9ny (7.2) 
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Only the exponential growth rate is significant here. We then apply Proposition [T2] to bound 
the number of zeros by the maximum of 

log I'lf ''f < CA + 0(logA) +0(log )■ 

\\OtUx\\L2(dn\\ \\OtUx\\L2{dn\\ 

Of course this estimate assumes u\{t) ^ identically, as can happen with radial eigenfunc- 
tions on the disc. Assuming (11. 2p . the third term of (17. 3p is 0(A), completing the proof. 
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